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I. INTRODUCTION 

The overall objective of this grant is to further the industrial 
applications potential of ion-beam etching and deposition processes. 

To accomplish this overall objective, a variety of specific tasks have 
•been undertaken. This report is divided into sections corresponding to 
the various tasks. 

Extensive experimental studies of the ion beam texturing of sur- 
faces were begun during the last support period^ and have been con- 
tinued and expanded. Primary emphasis has moved toward experimental 
checks on a surface texturing model developed during this past year. 

An end objective of this task has been to develop convenient and 
reliable procedures for texturing a variety of surfaces. The textured 
surfaces could be used for functions such as the control of secondary 
electron emission, the Increased absorption of light, and the Improved 
attachment of biological tissue (for medical Implants). 

Textured surfaces, typically with conical structures, have been 
produced by simultaneously etching a surface and seeding that surface 
with another material. A model based on surface diffusion predicts a 
variation in cone spacing with surface temperature, as well as a criti- 
cal temperature below which cones will not form. Substantial agreement 
with the model has been found for several combinations of seed and sur- 
face materials. Including one with a high sputter yield seed on a low 
sputter yield surface (gold on aluminum). Coning with this last com- 
bination was predicted by the theory for a sufficiently mobile seed 
materials The existence of a minimum temperature for the foirmation of 
cones should also be Important to those Interested in ion-beam 
machining of smooth surfaces. Elements contained in environmental 


contaminants or within the sputtered alloys or compounds may also serve 
as seed material. 

A mathematical model has been developed to calculate the sputter 
deposition of seed material onto a substrate during texturing. The 
deposition rates calculated using this technique are used as Inputs to 
the predictive model for cone formation and distribution. 

Calculations have been continued on cross sections for energy and 
momentum transfer in binary collisions between like pairs of Ar, Kr, 
and Xe atoms in the energy range from about leV to lOOOeV. These cal- 
culations were begun during the last support period using a classical 
model with pair interaction potentials for the rare gases developed 

from experimental sources, e.g., investigations of specific heats, 

2 

viscosities, solid state parameters and scattering data. Cross sections 
in this energy range have not been available. The cross sections 
exhibit a rapid decrease from accepted values at thermal energies as 
the interaction energy increases. This behavior can be used to under- 
stand directed beam propagation limits as well as particle propagation 
within a plasma. 

Silicon has been successfully textured on a micron scale with a 
higher activation energy being indicated than is typical of metallic 
systems. Operation of an ion source and sputter etch rates have both 
been studied using reactive gases instead of Ar, the usual working gas. 
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II. ION BEAM TEXTURING OF SURFACES 

Ion beam sputtering is frequently used in surface etch or milling 
operations where the desired result is the uniform removal of material, 
perhaps using a mask to produce a relief pattern on a substrate, or 
simply as a thinning process such as sample preparation for transmission 
electron microscopy. The processed surfaces, however, are often not 
microscopically smooth but exhibit a variety of surface textures that 
develop as a result of the ion bombardment. 

Tgxtures can develop on solid surfaces in several ways during 
directed ion beam sputtering. Physical sputtering of solid materials, 
which is customarily accomplished using chemically inert ion species, 
f>an reveal individual crystal grains because of differential sputter 
etch rates dependent on crystal grain orientation. Impurities, pre- 
cipitates, phases or initially irregular surfaces can also cause char- 
acteristic textures to develop during sputtering because of preferential 
of higher sputter yield sites. In addition, reactive ion beam 

etching can chemically texture a surface through selective removal of 

2 

elements from an alloy. 

Of particular interest here, however, is the texturing induced by 

the deliberate deposition of an Impurity onto a solid surface while 

simultaneously bombarding the surface with an ion beam. This technique 

is often referred to as '^seeding” , with the impurity being termed the 

seed material. Under appropriate conditions, microscopic cones or 

hillocks develop because of preferential sputtering of surrounding 

material. It has been generally understood that these cones result 

from clusters of seed atoms protecting the underlying substrate while 

3 

surrounding substrate material is etched away. 
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Ion beam texturing has been attempted with many material combi- 
nations and has several potential applications. Textured surfaces have 

been successfully used for enhanced absorption of radiant energy in 

4 

solar collectors. Perhaps one of the more promising applications is 

in the realm of biomedical materials such as pros theses with soft or 

5 6 7 8 

hard tissue interfaces that require firm bonding. ’ 

It has been postulated that the formation and replenishment of 
seed clusters is a result of the surface diffusion or migration of seed 
atoms with nucleatlon or attachment occurring when other seed atoms are 
encountered. It has also been widely believed that a necessary, but 

possibly not a sufficient, condition for the formation of sputter cones 
is that the seed material must have a lower sputter yield than the sub- 
strate material to account for the observed differential etch. The 
validity of these ideas is explored within the context of an analytical 
model along with experimental tests. 

A model has been developed as a first approximation to both a 
qualitative and quantitative description of the surface migration of 
seed atoms, the clustering process, and the resulting distribution of 
cones on a substrate. Experimentally controllable parameters as well 
as physical properties of the materials involved are included in the 
model. The diffusion aspect is described first. The critical minimum 
size of seed clusters is then related to the seed diffusion rate. As a 
further result, a critical substrate temperature is obtained, below 
which cones will not be formed by the diffusion and sputtering of seed 
material. The predictions based on this model are investigated experi- 
mentally and an assessment la made of the Impact of the model on the 
variously held postulates and assumptions regarding the seeding and cone 
formation processes. 


Sputter Cone Seeding Theory 

The primary features of a conceptual model used to understand seed 
diffusion are: 1) a substrate surface characterized as an array of 

uniformly distributed adsorption sites at which seed atoms can be bound 
to the surface after being deposited onto the substrate; 2) these 
indiyidual adsorption sites described as potential energy minima that 
interact with seed atoms through isotropic binding and restoring forces; 
3) the seed atoms of sizes and masses appropriate to the seeding species 
used that are deposited onto the substrate through some external means; 
and 4) an ion beam of a given energy and type of ion providing a uniform 
flux density of energetic ions over the substrate. 

Surface Diffusion 

The frequency for an isolated seed atom, located at an adsorption 
site and belonging to a population of adsorbed seed atoms in thermal 
equilibrium with a substrate of temperature T, to jump between 
adsorption sites is 

V = C1/Tq) exp (-Ej/kT) , (2-1) 

where is the characteristic time for the jumping process, is 

the activation energy for diffusion (i.e. , the energy barrier between 
adsorption sites) and k is Boltzmann's constant. The exponential 
factor reflects the probability that, given thermal equilibrium at tem- 
perature T, a seed atom has an energy greater than E^, assuming a 
uniform density of states. The reciprocal of the characteristic time 
may be thought of as a characteristic frequency of oscillation for 
the seed atom in the potential well representing the adsorption site. 


The total number of adsorption sites visited by an adsorbed seed atom 
during the time that it spends on the surface is given by the Jump fre- 
quency V multiplied by the mean adsorption time t : 

3i 

t>a = Tg = (Tg/TQ) exp (-Ejj/kT) . (2-2) 

The magnitude of the mean adsorption time is governed by those 

proee«f^es that would tend to remove or desorb a seed atom from the sur- 
face. One candidate process would be evaporation or sublimation which 
would occur if a seed atom possessed simultaneously sufficient energy 
to break the adsorption bond and a momentum directed away from the sur- 
face. If the seed atoms are considered to be in thermal equilibrium 
with the substrate, the probability of having sufficient energy to 
desorb would be exp(-E^/kT) where E^ is the total binding energy 
of the seed atom to the surface. Clearly, Eq>Ej, making this process 
highly Improbable. A more likely mechanism for removal of adsorbed 
seed atoms is physical sputtering by an applied high current density 
ion beam. The probability of a seed atom being sputtered depends 
linearly on the ion flux and on the seed cross section. 

The migration of seed atoms along the surface of the substrate can 
be viewed as a random-walk process. Using a^ as the mean distance 
between adsorption sites, the random-walk diffusion length r^ is 
given by 

’'d = Va*^ = (-Ed/2kT) . (2-3) 

The random-walk diffusion length represents the average net separation 
distance of a seed atom from its Initial adsorption site, achieved 


during a time t . This, then. Is the average radius through which 

seed atoms are expected to be able to diffuse. This expression (Eq. 

(2-3)) Is similar to diffusion equations derived for condensation of si 

vapor phase onto a substrate. It Is necessary to obtain appropriate 

values for x and x to apply this equation to the seeding and tex- 
a o 

turlng problem. 

The seed adsorption time Is limited by the Ion beam sputter 
removal of seed atoms once they are attached to the surface. If the 
Ion flux (ions/m^-sec) Is given by and the cross section that a 

seed atom presents to the Ion beam Is given by a then the probability 

SL 

per unit time P for Individual Ion-seed collisions to occur Is 
c 

- R.o^ . (2-A) 

c la 

The probability that a given collision will result in removal of an 
atom Is the sputter yield or sputter probability Y . Thus, the proba- 

8i 

bility per unit time P for sputter removal of individual seed atoms 

s 

from the substrate Is 


P 

s 


Y R.o 
ala 


(2-5) 


The reciprocal of this probability per unit time Is then the mean 
adsorption lifetime for a seed atom on the substrate 



1/Y R.a 
'a 1 a 


( 2 - 6 ) 


A first estimate of x can be obtained by inserting approximate values 

a 

for the quantities in Eq. (2-6), For a loosely bonded and mobile seed 
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at™, the aputterlng yield should be near unity. The crons section 

a, can he approximated as N„-2/3. „Heee is the atomic density of 

the substrate with a typical value of 10«„-3.9 1^3 


,19 


a typica^O A/m2 (1 mA/cm^) Ion beam current density Is about 6 x 10 

(m -sec) ^ with these substitutions an adsorption time of about 1/3 
sec is obtained. 

seed atoms of interest In the process of sputter cone formation 
must be mobile. The seed atoms should therefore exhibit weaker bonding 
to the substrate than that of the substrate atoms to one another. 
Because of this relatively weak bonding for seed atoms, the character- 
istic time fox a seed atom to jump between adsorption sites should 

be a function of this bonding Interaction, rather than the substrate 
lattice vibrations to which it Is only loosely coupled. Considering a 
seed atom at ar adsorption site, the site can be thought of as a poten- 
tial well within Which the seed atom Is bound. The potential barrier 
between adjacent sites Is the activation energy E,. Motion along the 

surface should therefore be characterised by a well depth E^ and an 

approximate periodicity, bslng a„ as the average spacing between 

sites, a sinusoidal representation of a periodic, undulating potential 
is 


U(r) ^ Ej ( 1 -cos 2 irr/a )/2 . 


( 2 - 7 ) 


A parabolic potential results from retaining only the first two terms 
in the cosine series expansion. 


U(r) % Tr^E,r^/a ^ 
d o 


( 2 - 8 ) 


Using as the mass of the seed atom, the equation of motion is 


m d^r/dt^= -au(r)/3r = -27T^E,r/a ^ , 

s do 


(2-9) 


which Implies an harmonic oscillation frequency of 


V = (E,/2m a 

d so' 


(2-10) 


or a period of 


X = 


(2m a /E.) 
s d 


h 


( 2 - 11 ) 


This expression probably underestimates somewhat the time to jump 

between adsorption sites because of the parabolic approximation that 

was assumed. A numerical estimate of x^ can be obtained using the 

npproprlate substitutions. The distance a can be approximated as 
-1/3 

» with again the atomic density. The mass m^ of a repre- 

sentative heavy seed atom can be taken as 3 x 10~^^ kg , while the 
activation energy E^ may be assumed to be in the neighborhood of 

1.6 X 10 J (leV).^^ With these substitutions, x is found to be 

o 

-13 

about 4 X 10 sec. 

Using Eqsn (2-6) and (2-11) to substitute for x and x in 

a o 

Eq. (2-3); 


'd ■ (-E^/2kI) 


(2-12) 


With the substitutions of N and N for a and a 

o o o a’ 


a simplified expression is obtained, 


together with unity for Y 


(2-13) 


'd - /2M) 


The exponential factor varies rapidly over a vrtde range. The square- 
root factor is much more slowly varying, so that numerical values may 
be substituted therein for a first approximation. Making the substi- 
tutions of 1.6 X lo'l’ j for lO^V^ for 3 x lo'^^ ^ 

»3 and 6 X 10« („2-sec)-l for i„ the squarLoot factor, 


Tj = 2 X 10 
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^xp (-Ej/2kT) . (m) 


(2-14) 


This, then, is an estimate of the radius over which surface diffusion 
can be expected to take place. This radius is a function primarily of 
activation energy and the substrate temperature, other fac- 
tors such as seed atom mass m^, l„„ fi„r rhe adsorption 

Site spacing a^ are considerably less important. The activation 
energy is a function of the seed-substrate material combination 
(together with a small temperature effect) and ranges from about 0.5 to 
2eV for metallic materials of interest for seeding. The large varia- 
tion of tj with temperature is shown in Fig. 2-1 for this range of 
activation energies. It can be seen that relatively modest temperature 
changes can have a profound effect on the mean seed diffusion radius. 

A distance equal to twice the average diffusion radius can be used 
as a measure of the average separation between clusters. If clusters 
were, on the average, much farther apart than 2r^, more clusters 
would begin to nucleate and grow in the intervening spaces where seed 
densities were enhanced, thus narrowing the gaps between clusters. 


V 





Fig. 2-1. Diffusion radius as a function of activation energy and 
substrate temperature for a specific case. 


At the other extreme, if clusters began to nucleate at separations much 
less than 2r^, larger clusters would grow faster and intercept dif- 
fusing seed material at the expense of smaller clusters. Thus, an 
average cluster separation of about 2r^ would be expected to be 
stable. 

Seed Clustering 

Having foujid the radius from which diffusion will supply a seed 
cluster, it is appropriate to consider the stability of that cluster. 
There is a critical size of a seed cluster, below which steady growth 
is not possible. This critical radius is essentially the same value 
as is used in nucleation theory. It is obtained by setting to zero the 
derivative of free energy (surface plus volume) with respect to radius 
Continuous growth can take place above this radius, but dissociation of 
the cluster will tend to occur at smaller radii. 

Whether or not a seed cluster is stable will therefore depend on 
the difiiusion rate to the cluster being sufficient, or insufficient, to 
supply the sputtering loss from a cluster of critical radius. To 
determine the stability requirements of seed clusters it is necessary 
to investigate more closely the processes of seed movement, and clus- 
tering while adsorbed on a substrate. 

Seed atoms are assumed to move from adsorption site to adsorption 
site on the surface by a random walk process. Only those seed atoms 
that acquire an energy greater than are mobile. Energy is 

exchanged between seed atoms and the lattice through the loose coupling 
of the seed atom to thermal lattice vibrations bringing the adsorbed 
seed atom population into approximate thermal equilibrium with the 


substrate lattice. It is then the temperature of the substrate that 
governs the fraction of seed atoms that are mobile at any given time. 

A random walk process for seed motion requires that, on the average, 
seed atoms experience an inelastic "collision" as a result of each jump 
to a neighboring adsorption site. The energy loss in the inelastic 
process would then be associated with a loss of the initial direction 
of motion resulting in a random direction for the next jump. Regarding 
each jump to a neighboring adsorption site as an independent event ter- 
®ij^3tlng with directional and energy losses can be further supported by 
considering the adsorption bonding forces tending to accelerate a seed 
atom down into the neighboring adsorption site once it has cleared the 
potential barrier E^, where is the diffusion activation energy. 

It is improbable that the incoming sesfd velocity vector would be suf- 
ficiently well-aligned with the spatial symmetry of the neighboring 
adsorption site to proceed without significant deflection. 

To appropriately model the clustering of seed atoms, it is first 
necessary to characterize the mean speed for the surface migration of 
seed atoms. There is the possibility of assigning different character- 
istic velocities to mobile seed atoms. One velocity, designating the micro- 
velocity ^jjj> could be derived from the spacing of adsorption sites 
and the jump time. 



(2-15) 


or, using Eq. (2-11), 
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Substituting leV for and 3 x 10 Kg for nig yields, 

u = 5.2 X 10^ m/sec . (2-17) 

m 

Another velocity, designated macro*velocity can be obtained from 

the random walk diffusion radius and the mean adsorption time, 




or, using Eqs. (2-6) and (2-13), 


(2-18) 


“ (Ej/2mg)^(R^Yg^/N^)*^exp(-E^/2kT) . (2-19) 

19 -2 -1 

Using the representative values of leV for E,, 6 x 10 iri sec for 

d 

-2S 9Q 

R., 1 for Y , 3 X 10 Kg for m , lO^" m for N , and 800®K for 

X a s o 

T we get 


= 4*0 X 10 ^ m/sec • (2-20) 

This is nine orders of magnitude less than Perhaps a more appro- 

priate expression for the seed velocity would take into account the 
total distance of travel found by summing all jumps taken in a seed 

lifetime t . The number of jumps is given by Eq. (2-2) : 
a 

n = (Tg/x^) exp (-E^/kT) , (2-21) 


and the mean seed speed is: 
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'"s “ ’ (2-22) 

or, using Eqs. (2-21) and (2-11) 

Og = (Ed/2mg)*^ exp (-Ed/kT) , (2-23) 


or, 

-4 

Wg =2.6 X 10 ra/sec (2-24) 

with the usual substitutions. Equation (2-23) then, is felt to be the 
most appropriate expression for average seed velocity for the purpose 
of understanding diffusion and clustering. 

With an expression for seed velocity, the detailed movements of 
seed atom populations will now be considered. As a first example, an 
estimate for the upper bound of the areal density of seed atoms on the 
substrate, N^, can be obtained by assuming that clusters fail to 
form on a seeded surface. In the absence of seed clusters, a steady 
state will be reached with the seed deposition rate and the re— sputtering 
rete being equal. The seed flux (seed atoms deposited/m^— sec) can 

be written in terms of the ion flux using a factor F , the seeding 
fraction, that reflects the ratio of seed deposition to ion flux: 


Rs ^s^l • (2-25) 

The rate of seed removal by physical sputtering is the probability 

per unit time for sputtering an individual seed atom multiplied by the 


areal seed density 


= Y R.a N 
s alas 


(2-26) 


Equating the deposition and sputter removal rates gives, 


F R = Y R.a N 
SI alas 


max 


(2-27) 


where Is the maximum achievable areal seed density In the 

absence of clustering. Solving for R- : 

^max 

^®max ~ ^s^^a*^a * (2-28) 

With a value of 0.01 for F^, this expression yields, using the cus- 
tomary substitutions: 


Nc 


’max 


% 2 X 10^^ m"^ . 


(2-29) 


We can also look at the fraction of surface coverage by seed atoms. 
Since 1/a^ is approximately the density of adsorption sites we get 


F N 


s a 



(2-30) 


where Y^ has been assumed to be approximately 1 and N is the 

a 

adsorption site density. Thus, the maximum fractional surface coverage 

by seed material is about equal to the seeding fraction F • 

s 

A more detailed analysis is necessary when clusters are assumed 
to be present on the substrate. Figure 2^2 represents a cluster and its 
surrounding seed diffusion area. The various N’s shown in Fig. 2-2 
are flow rates of seed atoms in atoms/sec. Let be the areal density 
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of seed clusters, then 1/N^ = A = + A^, where A^ is the cluster 

area and A^ is the diffusion 'drawing' area. The definitions of the 
seed flow rates are in Table 2-1. 


Table 2—1. Seed Flow Rates 

Seed atom deposition rate into area A, 

N~: Seed atom deposition rate into area A 

^3 seed atom acquisition by the cluster 

Rate of re-sputtering of seed atoms from 
Nc Rate of re-sputtering of seed atoms from A 

j c 

Ng Rate of diffusion of seed atoms out of A^ 

N, Rate of diffusion of seed atoms into A, 

' d 


All rates, are positive with the direction of seed flow Indicated 

by the arrows in Fig. 2-2. 


A steady state situation will be considered first. Under steady 

state conditions, the following relations hold: 





• • 


= N3 


N2 + N3 = N3 . 


(2-31) 

(2-32) 

(2-33) 


(2-34) 



Equation (2-31) simply expresses the S 3 nranetry of seed diffusion into 

and out of the area A containing a single cluster. Equation (2-32) 

equates the total seed deposition rate to the total sputter removal 

rate in area A. Equation (2—33) relates the net seed flux into A, 

d 

(taking re-sputtering into account) to the diffusion flux from A, to 

d 

the cluster. Finally, Eq. (2-34) equates the total flux to the cluster 
to the total sputtering from the cluster. Equation (2-32) can be 
obtained by eliminating between Eqs. (2-33) and (2-34). 


Most of the individual rates can be obtained from consideration of 


simple sputtering rates : 



R.F A, 
X s d 


(2-35) 



R.F A 


1 s c 


(2-36) 


N, = R.Y a N A^ 
4 X a a s d 


(2-37) 



R.Y A 
i c c 


(2-38) 


where Y^ is the sputter yield of atoms from the seed cluster. is 

related to the Vcollision’ frequency of seeds with clusters. Assuming 
that each seed colliding with a cluster adheres we can use a length of 
2rc to describe the two-dimensional cluster interaction cross section 
where r^ is the cluster radius. 


A 


- 20 - 


Thls relation can be understood more simply by noting that 2r u N 

c s c 

Is the collision frequency for a single seed moving through clusters of 
density at a speed u^. Multiplying by the number of seed atoms 

per unit area gives the total frequency of seed-cluster collisions per 
unit area. Multiplying by the "drawing" area for a single cluster 

gives the number of collisions with one cluster per unit time. 

With substitutions for the various rates, Eqns. (2-33) and (2-34) 
can be rewritten as: 


R^F A, -R Y o N A. = 2r u N N A, 
Isa laasd cscsd 


(2-40) 


and 


R.F^A^ + 2r u N N A, = R^Y A 
isc cscsd i c c 


(2-41) 


These equations may be solved simultaneously for the drawing area A^ 
and for the surface seed density N . 


A^(Y -F^)(2r u +A Y a R J 
A = ^ c s c s c a a 1 

d 2rcUgFg - A^(Y„-F„)Y„a„R, 


a a 1 


(2-42) 



A„(Y -F )R. 
c c s 1 

2r u 
c s 



(2-43) 


Under typical seeding and sputtering conditions the additive terms in 
Eq. (2-42) involving the ion arrival rate can be neglected giving: 


(2-44) 


Physically, this approximation corresponds to the actual seed coverage 
between clusters being much smaller than the maximum possible coverage. 
Equation (2-44) can be written in terms of r^ assuming circular 
geometry. 

'd ■ (2-«) 

This is the required diffusion radius to sustain a seed cluster of 
critical radius r^. For most applications, so that may 

often be neglected in the numerator with little error. 

Critical Temperature 

This diffusion radius (Eq. (2-45)) can be equated to that of 
Eq. (2-12) to determine the minimum substrate temperature that will 
give the required diffusion to sustain clusters of the minimum size. 

The logarithmic factor will vary slowly with changes in the variables 

therein; thus, the activation energy will most strongly dominate the 

critical temperature. The new variables that have been Introduced are 

the critical radius r^ and the seeding fraction F . The seeding 

fraction is calculated in Section III, a value of 0.01 will be used as 

a reasonable approximation. Values of the critical radius range from 

about 6 to 10 A cross check was made on these values of r from 

c 

another point of view. The use of a critical radius derived from 
nucleation studies Implies perturbations corresponding to thermal 
energy. The energy of an ion impact in a sputtering environment, 
however, is far above the thermal level. The limiting condition for 


Ion impact is that the cluster should be large enough to absorb all the 
ion energy without vaporizing. The heat of sublimation for typical 
seed materials ranges from about 3 to 9eV/atom.^^ For typical ion 
energies of 500 to lOOOeV, then, this minimum cluster radius is about 
5 A to 10 A. The heat of sublimation is, of course, a macroscopic 
quantity. However, this calculation indicates that the increase in 
cluster size to accomodate the maximum perturbation of an ion collision 
should be moderate, almost certainly less than a factor of two times 
the nucleatlon theory radius of 10 A. 

Using 10 A for r^ and 0.01 for an approximate expression 

can be obtained for T^ in terms of the activation energy measured in 
eV .for this expression only. 


T^ Jb 587 “K E, 
c d 


(2-47) 


This expression may be used for rapid estimates and approximations 
without paying undue attention to the precise values used for some of 
the other quantities. 

Experimental Results and Discussion 

Two physical effects following from the cone seeding theory were 
investigated experimentally; ( 1 ) the variation in cone density with 
substrate temperature, and (2) the existence of a critical temperature 
for cone formation. In the experiments, samples of aluminum were 
heated to a specified temperature and sputtered for 20 min. using a 
1 mA/cm beam of 500eV Ar ions, while simultaneously seeding the sample 
with a source of either Mo or Au. The specified substrate temperatures 


ranged from 400 to 575®C in different tests, with In used to assure 
good thermal contact between the samples and the heat source. 

Seeding was accomplished by partially covering the sample with a 
sheet of the seed material. Seed material was then sputtered from the 
beveled edge of this covering sheet onto the sample. Because of this 
geometrical arrangement, the seeding intensity decreased with increasing 
distance from the beveled edge (see Section III) . After sputtering, 
scanning electron micrographs were obtained of the sample surfaces. 

The average separation between cones was then (A/N)^, where N is the 
number of cones counted in an area A of the micrograph. 

We identify <r>, one half the average separation between cones, 
with the diffusion radius r^ in the preceding theory. Figures 2-3 
and 2-4 show plots of <r> as a function of reciprocal temperature for 
the two seed materials. From the theory, we expect straight lines on 
these plots. This can be shown by rewriting Eq. (2-13) as 

= r^exp(-E^/2kT) (2-48) 

where r is t constant 
o 

■ ( 2 - 49 ) 

Taking the natural logarithm of Eq. (2-49) and rearranging the result, 

An r^ = (-E^/2k) 1/T + In r^ . (2-50) 

This expression is linear if In r^ is plotted as a function of inverse 


temperature, as done in Figs, 2-3 and 2-4. The slope of the line is 
-Ej/2k, while the ordinate intercept Is r^. 

A least squares regression line was fitted to each of Figs. 2-3 
and 2-4. The best fit corresponded to an E. of 1.20eV and an r 

a O 

7 o 

of 1.10 X 10 A for Fig. 2-3 and an E. of 1.04eV and an r of 

u O 

6 

6.04 X 10 A for Fig. 2-4. The correlation coefficients exceeded 0.9 
for both fits. Experimental critical temperatures were also determined: 
700 K for A1 seeded with Mo and 675°K for A1 seeded with Au. 

The qualitative trends shown in Figs. 2-3 and 2-4 are all con- 
sistent with theory. The data shown also include varying distances 
from the edge of covering sheet, which corresponded to different seed 
ratios F^. The only data selection used was that regions far enough 
from the covering sheet (low enough seed ratio) showed large fractions 
of the areas without cones and were omitted. For regions with sufficient 
seed to show fully developed cone structures, then, the effect of seed 
fraction on cone spacing must be small. This result is in agree- 

ment with Eq. (2-13), which does not include F^ as a variable. Acti- 
vation energies for surface diffusion were not available from literature 
for the exact combinations of materials used, but the values obtained 
from the slopes of Figs. 2-3 and 2-4 are in qualitative agreement with 
other combinations that have been studied using standard diffusion 
measurement techniques. 

The activation energies from Figs. 2-3 and 2-4 can be used for 
a more quantitative comparison of experimental results with theory. 

Using these activation energies together with 6.24 x 10^^ ions/m^-sec 
for *01 for F^, 1 x 10 m for r^, and handbook values for 

Other variables, theoretical values were calculated for r and the 


Mean Diffusion Radius, <r>, 
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1.15 1.20 1.25 1.30 1.35 

inverse Temperature, I /T, 


Fig. 2-3. Correlation for A1 seeded with Mo. 



Mean Diffusion Radius,<r> , A 



Fig. 2~4* Corrialatlort for A1 seeded with Au. 
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critical temperature. These values were 2.21 x 10^ A and 684®K for 
A1 seeded with Mo and 1.78 x 10^ A and 646®K for A1 seeded with Au. 

The agreement for critical temperature is close, with experimental 
values being 2 and 4 percent above the theoretical values from Eq. 
(2-46). This close agreement is due in part to most of the variables 
in Eq. (2-46) being included in a logarithmic factor. Closer agreement 
is to be expected between the predictions of the model and experimental 
results for the critical temperature than for other aspects of the 
model. For example, a limiting case of a flat surface as used in this 
model is more appropriate for derivation of the critical temperature 
at which no cones are present than for determining cone separation on 
a fully developed surface. The agreement for cone spacing is less 
close, with experimental values of r^ and cone spacing being 3 to 5 
times the theoretical ones from Eq. (2-13). This difference may be 
the result of theory being derived for a flat isotropic surface, while 
the experimental values were for well developed conical surfaces. The 
cascades of sputtered substrate atoms on the conical surfaces would be 
expected to reduce seed diffusion toward the apexes, resulting in the 
reduction of some seed clusters below the critical size. 

The choice of materials for Fig. 2-4 also represents support for 
the theory. Previous descriptions of coning states that low sputter 
yield materials protected the apexes of cones in higher sputter yield 
substrates. The present theory indicates that the seed can have a 
higher sputter yield than the substrate if the seed material is “,'f- 
fliciently mobile to replenish the sputtering loss from a seed cluster. 
The Au seed in Fig. 2-4 has a much higher yield than the A1 substrate, 
so that this result is an example of coning that was predicted by the 
present theory. 




A texturing theory based on surface diffusion has shown substantial 
agreement with experimental results. The existence of a minimum sub- 
strate temperature for texturing, the variation of cone spacing with 
temperatures and the texturing with a high sputter yield seed on a lower 
sputter yield substrate were all predicted by this theory. This theory 
has obvious utility in the production of textured surfaces for applica- 
tions such as low reflectance solar cells, high emissivity heat radia- 
tors, and medical Implants. The existence of a minimum critical tem- 
perature for coning is also ir.,portant for the production of smooth 
surfaces, in that a sufficient reduction in surface temperature should 
reduce the mobility of any seed material enough to avoid coning, 
whether the seeding is due to stray sputtering or one of the elements 
in the substrate material. Due to the continual cleaning of the surface 
by the. bombarding specie, the coning process also permits evaluation of 

the surface diffusion activation energy under conditions of extreme 
cleanliness. 
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III. SEEDING RATE PROM A BEVELED SOURCE 


Experimental arrangements for ion beam texturing include a beveled 
source of seed material placed on top of a substrate. Part of the sub- 
strate is covered with the seed source providing a control surface for 
comparison with the part that is uncovered and exposed to both sputter 
deposition of seed material and ion bombardment. Such an arrangement 
is shown in Fig. 3—1. This study will focus on the deposition rate of 
seed material onto the substrate from such a beveled surface. 

A number of assumptions will be implicit in this seeding rate 
model. The ion beam will be assumed to have uniform current density 
over the area of interest and the beam will be directed normal to the 
substrate surface. The emission of sputtered seed atoms from points 
on the beveled edge will be assumed to correspond to a cosine distribu- 
tion about the local surface normal. This assumption is in close agree- 
ment with experimental observations where only small departures from 
a cosine distribution have been observed. 

To model in detail the emission of sputtered atoms from the 
beveled surface and their subsequent deposition onto the substrate, we 
explore the geometry indicated in Fig. 3-2. Let dF 2 be the rate of 
sputtered particle collection at the element of area da 2 in particles 
per unit time arriving from the element of area da^. This rate is 
proportional to the fraction of hemispherical solid angle subtended by 
da 2 as viewed from da^^, i.e., da2COS(}i2/2TTr^. The emission from 

‘^^l ‘f’l is proportional to da^^coscj)^ because of the 

assumed cosine emission distribution. The emission from arriving 

at da 2 


Can thus be expressed as: 
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Fig. 3-2. Geometry for sputter emission and flux collection. 
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• , 


dV2 * F da^cos(|>^ 


d&2^os^2 


2Trr‘‘ 


( 3 - 1 ) 


where F is a constant of proportionality representing the sputtered 
flux in particles per unit time per unit area emitted normal to da^^. 
The factor F normalizes to total emission from the emitter plane, 

Integrating over the plane P, will yield the total rate of par- 

cos4>2 

— — cos(|»j^dai ‘ (3-2) 

r 

An effective flux at da 2 is found by dividing the particle arrival 
rate by the area of the element da 2 


tide collection at da^ 


Fda, 


dr« = 


2 2ir 


/ 


F 


2 


p r cos <|)2 

^ " 2ir 7 T~ cos(|)j^daj^ . 

"" 


( 3 - 3 ) 


The relative seeding rate as a function of position on the sub- 


strate plane can thus be found from the Integral indicated in Eq. (3-3). 
This integral will, in general, not be solvable in closed form and will 
thus require a numerical precedure for evaluation. 

Figure 3-3 shows the geometry and definition of variables to be 
used in a numerical Integration algorithm. From Fig. 3-3 it can be 
seen that: 

= (x2-‘Xj^)^ + y^^sln^3+ (y2+y2Cos6)^ , (3-4) 



Fig. 3-3. 


Geometry and parameter definition for numerical integration algorithm. 
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co 8*2 » yj^slnS/r , 


(3-5) 


and that 


cos((>j^ » y2sin0/r 


(3-6) 


Thus Eq. (3-3) becomes: 


P 


L h/sln3 2 

, F r ^ r ynyoSln 6 

2 ■ s ; '•=‘1 j «yi -H — • 


(3-7) 


Given Che coordinates of the observation point and the dimensions and 
bevel angle of the seed source Eq. (3-7) can be written! 


Fy2S±n^3 ^ h/sin3 


’'2 ■ 2r f !iyiyi[(!C 2 -ri)' 


2 2 " 
+ yj^ sin 3+ (y 2 + y^cos3)^ 


-2 


(3-8) 


The numerical integration is carried out in N increments 


M increments of where: 


of Ax, and 


AXj_ « L/N 


(3-9) 


and 


Ayj^ = h/Msln3 . 


(3-10) 
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The Integral can now be written approximately 


as a summation: 


where 


I2 

F 


y„sln e N M 

E E 

1=1 j=i 


h\) 


yn 




+ (yj+yijCosBjn 


(3-11) 


Xj = (i-%)Ax, 
1 ^ 


(3-12) 


and 


Fi - (3-^) Ay, . 
i + 


(3-13) 


A program was written to carry out the numerical procedure Indi- 
cated in Eq. (3-11). Ihe factor F la carried aa a constant through 
the numerical wort. It la, however, proportional to the ion arrival 
rate and the sputter yield of the seed material. 


F = cos3 Y(S) 


(3-14) 


where E^ la the Ion arrival rate In Iona per unit area per unit time 
at the substrate, the cosine of 6 takes Into account the tilt of the 
beveled edge with respect to the ion beam direction, and Y(B) Is the 
sputter yield of the seed material at the angle 6 In atoms per ion. 


First, the integration step size necessary to assure convergence 
within one percent was Investigated. The parameters used corresponded 
to a realistic experimental arrangement: L/h = 8.0, 8 = 45“ , X 2 /h = 4.0, 

and 72/h = The results of calculating F 2 /F as a function of N 

and M are given in Table 3-1. 


Table 3-1. Integrations with Various 
Increments of Area 


N 

M 

F 2 /F 


6 

1 

3.0645 

X 

10"2 

11 

2 

2.8013 

X 

10"2 

17 

3 

2.6867 

X 

10" 2 

23 

4 

2.6396 

X 

10"^ 

28 

5 

2.6164 

X 

10-2 

-2 

34 

6 

2.6035 

X 

10 

-2 

40 

7 

2.5955 

X 

10 

-2 

45 

8 

2.5903 

X 

10 

-2 

51 

9 

2.5867 

X 

10 

-2 

154 

27 

2.5746 

X 

10 

-2 

20 

^ 7 

2.5957 

X 

10 

-2 

10 

7 

2.5948 

X 

10 

-2 

8 

7 

2.5881 

X 

10 

-2 

6 

7 

2.5489 

X 

10 

-2 

4 

7 

2.3271 

X 

10 


Letting the N=154, M=27 value be the standard, this standard is 
approached within 1% by letting N = 40 and M=7. However, N can be 
further reduced by a factor of 4 without sacrificing significant 
accuracy. These values, N - 10 and M = 7, are appropriate when the 
observation point 72 /^ Ts a unit distance, or greater, from the source. 
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If Y 2 /h is decreased by a factor of two, accuracy can be maintained 
by increasing both M and N each by a factor of two. Other, similar 
adjustments would be appropriate as y 2 /h changes. 

The seeding rate variation as a function of distance from the 
center of the seed source is given in Table 3-2 for the choice of 
parameters: L/h =8.0, 3 = ^5®, X 2 /h = 4.0. 

The data of Table 3-2 are plotted in Figs. ■'3-4 and 3-5. The seeding 
rate decreases rapidly near the source and falls off more slowly at 
greater distances as shown in Fig. 3-4. For typical dimensions of seed 
source and substrate, the seeding rate can vary by one or two orders of 
magnitude over the substrate. Both seeding rate and distance are 
represented on logarithmic scales in Fig. 3-5. The linear portion at 
large distances from the source has a slope corresponding to an inverse 
cubic variation with distance in agreement with the J 2 dependence of 
Eq. (3-8). At small distances from the source, the seeding rate 
approaches a limiting value that would correspond to the seeding rate 
from a source of semi-infinite extent tilted at an angle 3« 

The limit can be calculated by integrating Eq. (3-8) after changing 
to suitable dimensionless variables: u = Xj^/y 2 and u = y 2 ^/y 2 * 

00 00 

-00 o ■ ■ ■ 

+ (14- ucos3)^J (3-15) 

Equation (3-15) can be Integrated to yield: 

(F2/F)j^imit = a-cos3)/4 (3-16) 




iT 



-38- 


Table 3-2. Seeding Rate Along the Centerline of the Source. 


72/11 ^2^^ ^2^^ 


0.03125 

7.07 

X 

10“^ 

2.6 

i 8.02 

X 

10“3 

0.0625 

6.81 

X 

10"^ 

2.7 

7.55 

X 

10-3 

0.125 

6.32 

X 

10-2 

2,8 

7.11 

X 

10-3 

0.5 

4.18 

X 

10"2 

2.9 

6.71 

X 

10-3 

0.6 

3.77 

X 

10“2 

3.0 

6.33 

X 

10-3 

0.7 

3,41 

X 

10"2 

3.1 

5.98 

X 

10"3 

0.8 

3.10 

X 

10"2 

3.2 

5.66 

X 

10-3 

0.9 

2.82 

X 

10"2 

3.4 

5.07 

X 

10-3 

1.0 

2,59 

X 

10"2 

3.5 

4.81 

X 

10-3 

1.1 

2.38 

X 

10"2 

3.6 

4.57 

X 

10-3 

1.2 

2a8 

X 

10"2 

3.7 

4.34 

X 

10"3 

1.3 

2,01 

X 

10"2 

3.8 

4.12 

X 

10-3 

1.4 

1.85 

X 

10"2 

3.9 

3.92 

X 

10-3 

1.5 

1.71 

X 

10~2 

4.0 

3.73 

X 

10-3 

1.6 

1,58 

X 

10“2 

8.0 

8.70 

X 

10“^ 

1.7 

1.47 

X 

10"2 

16 

1.37 

X 

IQ-'^ 

1.8 

1.37 

X 

10-2 

32 

1.83 

X 

10"3 

1.9 

1.27 

X 

10-2 

64 

2.35 

X 

10-^ 

2.0 

1.19 

X 

10-2 

128 

2.99 

X 

10-' 

2.1 

1.11 

X 

10-2 

256 

3.77 

X 

10-« 

2.2 

1.04 

X 

10-2 

512 

4.72 

X 

10-^ 

2.3 

9.70 

X 

10-^ 

1024 

5.92 

X 

0 

1 

0 

r— 1 

2.4 

9.09 

X 

10-3 

2048 

7.40 

X 

10-^3- 

2.5 

8.54 

X 

10-3 

4096 

9.26 

X 

10-^2 
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Fig. 3-5. Seeding rate versus distance from the seed source along the 


centerline (log-log plot). 




which can be evaluated for 3 = 45® as 7.32 x lO"^. This value is 
plotted in Fig. 3-5 as the y 2 /h = 0 limit which the function 
approaches very closely at small y 2 /h. 

Table 3-3 lists calculated values for F 2 /F obtained over a 
regular mesh covering the substrate region for the same geometrical 
parameters used for Table 3-2. The data of Table 3-3 are plotted in 
Fig. 3-6. Close to the seed source there is a strong variation with 
lateral position dropping to only about one half the center value at 
the extreme ends where flux contributions are being received from the 
right or left side only. As the observation line is moved further 
from the source the lateral variations diminish considerably because 
the subtended angle of the seed source is then changed less when 
observed from various lateral poaitlons. 

An Auger surface composition analysis done at various distances 
from the shield location on a seeded sample should provide an experi- 
mental verification of the calculated profile shown in Fig. 3-4. 

From seeding theory (see Section II) we have shown that the 
critical temperature varies inversely as the logarithm of the seeding 
density times a constant. 

11600®K 

"^c InF ^ ^ ^ ^ (3-17) 

' S.' 

for reasonable magnitudes of the variables (E^ in eV) . Table 3-4 shows 

a possible variation of T with F . 

c s . 
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Table 3-3. Calculated Values for F 2 /F for Positions on the Substrate. 


X2/h 

0 . 

5 

0.0 

2.11 

10 "^ 

0.5 

3.32 

10 "^ 

1.0 

3.82 

10-2 

1.5 

4.02 

10"2 

2.0 

4.11 

CM 

1 

0 

rH 

2.5 

4.15 

10”2 

3.0 

4.17 

10"2 

3.5 

4.18 

10"2 

4.0 

4.18 

0 

1 

N > 


73/^ 

1.0 1 


1.33 

10“2 

8.92 

1.86 

10"2 

1.16 

2.20 

10~2 

1.37 

2.39 

10“2 

1.51 

2.49 

10"2 

1.60 

2.55 

10-2 

1. 66 

2.58 

10"2 

1.69 

2.59 

10"2 

1.71 

2.59 

CM 

1 

0 

1.71 


•5 2.0 


10"2 

6.35 

10-2 

10“2 

7.89 

10-2 

10"2 

9.20 

10-2 

CM 

1 

0 

rH 

1.02 

10-2 

10“2 

1.09 

10"2 

10"2 

1.14 

10“2 

0 

1 

1.17 

10"2 

10"2 

1.18 

10-2 

10-2 

1.19 

10-2 


X2/h 


2.5 

0.0 

4.72 

10-2 

0.5 

5.67 

10-2 

1.0 

6.53 

10-2 

1.5 

7.22 

10-2 

2.0 

7.75 

10"2 

2.5 

8.12 

10-2 

3.0 

8.36 

0 

1 

U ) 

3,5 

8.49 

10-2 

4,0 

8.54 

c *^ 

1 

0 


72 /h 
3.0 


3.62 

10-2 

2.85 

4.25 

0 

1 

U) 

3.28 

4.83 

10-2 

3.69 

5.32 

10-2 

4.04 

5.70 

10-2 

4.32 

5.99 

10-2 

4.54 

6.18 

10-2 

4.69 

6.30 

0 

1 

4.78 

6.33 

0 

1 

4.81 


.5 


4.0 

10-2 

2.29 

10-2 

10-2 

2.60 

10-2 

10-2 

2.89 

10-2 

10-2 

3.14 

10-2 

10“2 

3.35 

10-2 

10-2 

3.52 

10-2 

10-2 

3.64 

10"2 

10-2 

3.71 

10-2 

10-2 

3.73 

10-2 



1.0 




Table 3-4. Variation of Critical Temperature 
with Seed Fraction. 


T 

c 

F 

s 

AT 

c 

735®K 

.0001 


641»K 

.001 

94“K 

569”K 

.01 

72“K 

511“K 

.1 

58°K 

464»K 

1.0 

47“K 

Thus, the variation in seeding density 

across a substrate can be trans 

lated into a corresponding 

variation in 

the critical texturing tempera- 

ture across the substrate. 

Contours of 

equal critical temperature can 

In theory be specified as to shape and position. If a substrate were 


held at a temperature corresponding to a contour a certain distance 
from the seed source, then locations between the source and that contou 
would be textured and locations beyond would remain non-textured. Sue 
a texture boundary would be convex, curved and symmetrical about the 
centerline of the seed source, silicon seeded with Mo has exhibited 

textured regions suggestive of this effect. 


The technique presented here for calculating seeding rate can be 

used to determine appropriate seed source geometries to be utilized 

when uniform seeding Is desired over an entire substrate to generate 
uniform texturing. 


IV. ENERGETIC BINARY COLLISIONS IN RARE GAS PLASMAS 


Energetic binary collisions become important in plasmas where 
high voltages are introduced to accelerate ions toward a target or 
to form a beam. Interest in these collision cross sections arose pri- 
marily in connection with experimental work involving the propagation 
of ion beams used for sputtering and processes occurring within diode 
plasma sputtering devices. Thus, although they should also be of 
interest in other areas, the cross sections obtained herein are 
largely described in terms of their Impact on the understanding of 
these particular experimental sputtering applications as representa- 
tive systems. Energies of several hundred eV are typically employed 
in diode sputtering and in ion beam devices.^ At such energies, 
interaction cross sections differ sharply from their values at very 
low or thermal energies.^ 

The plasmas and beams generated in sputtering devices are typically 
derived from rare gases with Ar by far being the most common. Other 
more reactive gases are employed occasionally to achieve special effects 
due to chemxcal reactions. This work is restricted to consideration of 
Ar, Kr, and Xe for momentum and energy transfer calculations. 

While it is recognized that electrons will generally be present in 
plasmas and propagating beams, their collisional effects will be 
neglected. Electrons are Important in plasmas and ion beams for space- 
charge neutralization and for current neutralization of an ion beam, 
but do not otherwise enter significantly into reactions with high 
energy ions. Electron masses are on the order of 10~^ of atom masses 
so that energy and momentum losses to electrons are negligible. The 
process of recombination of electrons with ions can likewise be 
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neglected teceuee the crose section for this process is several orders 
of magnitude smaller than other cross sections of Interest.’ 

Work on these interaction cross sections was begun during the 
previous support period.’ The work on resonance charge exchange will 
not be repeated here nor will the details of the calculation techniques 
The collision cross sections reported earlier however, dealt only with 
Ar and only with the cross section for directed momentum loss. Addi- 
tional cross sections for Kr and Xe have been calculated and two types 
of total cross sections have been addressed: momentum transfer and 

energy transfer, corresponding to the customary diffusion and viscosity 

cross sections respectively. 

Elas tic Collision Dynamics 

Only through scattering from particles of comparable mass can 
energetic particles lose both significant energy and forward momentum. 
This study Involves, specifically, elastic collisions between like 
neutral particles, but the qualitative trends may be useful lu under- 
standing ion-atom collisions or unlike atom collisions as well. The 
probability of charge exchange prior to elastic collisions is high for 
beam Ions. Thus, the collision of high velocity neutrals with low 
velocity neutrals is the dominant process for momentum and energy loss 
by beam particles prior to sputtering. 

The dcBroglle wavelength of an Ar. Kr, or Xe atom at typical ion 
energies Is very small compared to atomic dimensions, allowing the use 

of classical particle traiectorles at a f-tv-oi- a 

j ccorxes as a first approximation to the 

scattering. 

If the interaction potential between two atoms Is known, the clas- 
sical equations of motion can be solved to relate the scattering angle 


to the impact parameter for the Incident atom. The solution is most 

directly obtained using center of mass coordinates where a particle of 

reduced mass p interacts with a fixed center of force through a 

potential V(r) where r is the separation distance of the particle 

from the center of force. ^ Solving the equation of motion yields an 

expression for the scattering angle in the laboratory frame of ref- 
7,8 

erence 



where is the speed of the incident particle, s is the impact 

parameter, and r^ is the distance of closest approach. 

The integration of Eq. (4-1) is not straightforward because the 
upper limit is infinite and the integrand itself has a singularity at 
the point of closest approach. Closed form solutions in terms of 
elliptic integrals are known for only a few power-law functions of r 
for the potential V(r)-' The actual integration can be carried out 
numerically after a change of variables to u = 1/r which gives finite 
limits but does not remove the singularity in the integral. An algo- 
rithm was developed for this integration utilizing a ten-point Gauss- 
Legendre quadrature scheme along with a logarithmically decreasing 
integration interval as the integration proceeds from the well-behaved 
end of the interval toward the singularity. Before integrating, the 
upper limit of integration, was first determined by an 

iterative solution to the (usually) transcendental equation: 
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2V(r^)/E^ + s2/r^2 , ^ 


where Is the energy of the Incident particle In the laboratory 

frame of reference. 

To verify that the Integration algorithm was functioning properly, 
scattering angles for an Inverse square law force (Rutherford Scattering) 
were calculated. These calculated angles were compared with the closed- 
form solution and the numerical values exhibit a divergence of only a 
few units in the third significant figure. 

Interaction Potentials 

To use a classical approximation for potential scattering it is 
necessary to obtain a potential function describing the pair inter- 
action. Barker. Fisher and Watts^ have developed an Ar pair potential 
that is recognized as being in best overall agreement with experimental 
data from a number of sources such as specific heats, internal energies, 
vxscosities, solid state parameters, and scattering data. Barker et. 
al.^0 have developed potentials for Kr and Xe using nearly the same 

analytical form for the potential as for Ar with similar adjustable 
parameters. 

The screened Coulomb potential used by Bingham^^ and Everhart et. 
al. is valid at higher energies (e.g., above 50KeV) but disregards 
the electron repulsion important at lower energies. Smoothly joining 
the Bingham^ and Barker® potentials between the energies where each 
potential Is known to fit experlnental data well, supports the use of 
Barker's results up to about IKeV. 


Figure 4—1 compares the pair potentials used for Ar, Kr and Xe 
showing the increase in effective radius of the atom as the atomic 
number increases. The general character of the like-particle scattering 
from Kr and Xe in the classical approximation should be similar to Ar 
since the potential shapes are similar to Ar. However, the corres- 
ponding cross sections for the heavier gases should be larger at the 
same collision energy. 


Calculated Cross Sections 

Numerical solutions to Eq. (4-1) allow the calculation of differ- 
ential collision cross sections for different energies and gases. The 
differential cross section is given by 


a(e) = 


s ds 
sin0 d0 


(4-3) 


In a collision in which an atom is scattered through an angle 0, it 

1 ' ■ 

will have lost a fraction of its energy equal to (1-cos 0). Integratxng 

this fraction weighted by the differential cross section gives the cus- 

14 

tomary form for the widely used cross section for viscosity 




(l-cos^0) o(0) 2irsin0d0 


(4-4) 


A similar expression is obtained when the loss of total momentum is 
considered. This momentum loss cross section is used to describe dif- 
fusion processes. 




(1-COS0) o(0) 


2irsin0d0 


(4-5) 


Pair I 



Integration of Eqs. (4-4) and (4-5) yields the total cross sections for 
energy and momentum transfer given in Table 4-1 and plotted in Fig. 4-2 
and 4-3 over the energy range from leV to lOOOeV. 

Conclusion 

Rare gas interaction cross sections have been obtained that will 
allow a more thorough analysis and understanding of a variety of ener- 
getic processes involved in ion beam and plasma devices. Cross sections 
at these high energies have not previously been available to investi- 
gators, resulting in the occasional erroneous use of thermal values. 

The argon energy transfer cross section has been successfully applied 
to the elucidation, by Harper et. al. of a complex negative ion sput- 
tering phenomenon discovered in a diode plasma device. These cross 
sections for elastic collision processes along with the charge exchange 
cross section allow designers of ion beam deposition and machining 
experxments to quickly arrive at appropriate geometries and pressure 
environments for beam propagation. 



Table 4-1. Calculated Energy and Momentum Transfer Cross Sections. 


Energy, 

eV 

Energy Transfer 
Cross Section., 

Momentum Transfer 
Cross Section, A2 


Ar 

Kr 

Xe 

Ar 

Kr 

Xe 

1 

11.2 

14.0 

17.7 

7.17 

9.04 

11.5 

2 

9.54 

12.5 

15.9 

6.07 

8.10 

10.3 

3 

8.57 

11.8 

14.9 

5.42 

7.61 

9,66 

5 

7.37 

10.9 

13.8 

4.63 

7.03 

8.93 

10 

5.88 

V.78 

12.4 

3.68 

6.31 

8.02 

20 

4.72 

8.77 

11.1 

2.97 

5.66 

7.18 

30 

4.20 

8.22 

10.4 

2.65 

5.30 

6.73 

50 

3.67 

7.55 

9.59 

2.32 

4.86 

6.17 

100 

3.09 

6.70 

8.51 

1.96 

4.30 

5.47 

200 

2.60 

5.89 

7.49 

1.65 

3.77 

4.81 

300 

2.34 

5.44 

6.93 

1.48 

3.48 

4.44 

500 

2.03 

4.89 

6.24 

1.28 

3.12 

3.99 

1000 

1.65 

4.19 

5.36 

1.04 

2.66 

3.41 
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Fig. 4-2. Energy transfer (viscosity) cross section for Ar, Kr, and 
Xe as a function of energy. 
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3. Momentum transfer (diffusion) cross section for Ar, Kr, and 
Xe as a function of energy. 
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V. PHYSICAL SPUTTERING AND REACTIVE ION ETCHING 
by Curtis Haynes 

Experiments in seeding and surface texturing have been performed 

with an Ar ion beam using various combinations of seed material and 

substrates. In addition, the characteristics of ion source operation 

with CFy have been studied. 

4 

Ion Source Operation with CF^ 

A multipole ion source 7.6 cm in diameter and 5.1 cm long was 
tested using CF^ as the working gas. A refractory metal cathode was 
mounted in the center of the discharge chamber. A set of dished Mo 
grids was masked down to a beam diameter of 5 cm using 0.025 mm Ta 
foil. 

This source was operated with both CF^ and Ar for comparison of 
operating conditions. Figure 5-1 is a plot of discharge loss per 
extracted beam ion as a function of total beam current. The CF^ dis- 
charge tended to extinguish at low discharge voltages so it was operated 
at 75 V while Ar discharges are typically 50 V. Gomparlsons were made 
at different pressures with the pressure being measured at the back- 
ground pressure of the vacuum system and not inside the ion souroe 

-.4 ' ■ 

itself. Neither gas gave a stable discharge below about lx 10 Torr. 
From Fig. 5-1 it can be seen that the power requirements for CF^ are 
about 50 percent greater than for Ar for the same beam current* 

Looking at discharge current alone rather than discharge power 
(Fig* 5-2) it can be seen that the ratio of discharge Gurrcnt to beam 
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Pressure Ar 
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□ 

6 xlO’^Torr 
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0 

3xicT^ 
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A 

1 X 10’^ 


Btom Current, mA 

Fig. 5-1. Comparison of Ar and CF, source operating characteristics 








Fig* 5-2 i Correlation f 


or Ar and CF^ source operation. 




current correlates approximately linearly with beam current for both 
working gases. 

If the source is operated using CF^ with a constant cathode heater 
current there is a steady decrease in discharge current with source 
operation becoming impossible within about three hours. This reduction 
in discharge current is felt to be the result of a slow buildup of a 
non-conducting layer on the anodes, and a contributor to the data 
scatter in Fig. 5—1. The Ta cathode is also degraded by operation 
with CF^, becoming white and brittle within less than one hour of dis- 
charge operation. The 0.025 mm Ta foil used to mask the grids was 
eroded completely through with about ten hours total of CF^ discharge 
operation. In contrast, Ar used in the chamber for over 30 hours did 
not cause measurable erosion damage to the foil. 

Differences in the physical appearance of the interior of the dis- 
charge chamber were noticed after using each of three different pro- 
pellants. When operated with Hg the stainless steel surfaces developed 
a very shiny appearance that lasted through months :>f storage. In 
addition, the Fe pole pieces did not show any evidence of corrosion, 
whereas pole pieces that were made at the same time but never used 
inside a discharge chamber developed a layer of rust in the same storage 
environment. When operated with Ar, the Interior surfaces lost their 
shine and became a light, dull gray. Finally, a darker gray coating is 
observed after CF^ operation. 

Seeding and Texturing Experiments 

A variety of seeding and texturing experiments were done using a 
500eV Ar ion beam. Unless otherwise specified, the beam current 
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density was set at 1 mA/cm . Substrates were sputtered while resting on 
a heater with the temperature monitored using a thermocouple. 

Aluminum 

A1 was sputtered at 475°C using Fe as a seed material. Experiments 
have been performed using Fe as both a seed material and as a substrate 
in attempts to observe ferromagnetic clustering phenomena made possible 
by surface diffusion. The only system so far to exhibit what are 
believed to be ferromagnetic effects is an A1 substrate seeded with Fe. 
An array of long, narrow clusters of material could be seen with the 
long clusters preferentially aligned along what appear to be two 
mutually orthogonal directions (Fig. 5-3). It appears that the clusters 
might be ferromagnetic and aligned in preferred directions defined by 
the crystal planes of the underlying A1 which has an fee structure. It 
would seem reasonable that ferromagnetic coupling is responsible for 
the linear growth of the clusters because of the magnetic anisotropy 
that would be encountered by a diffusing Fe atom approaching a mag- 
netized cluster. 

In support of the hypothesis that these linear structures are 
aligned in some way with the underlying substrates, sharp boundaries 
could be located across which the direction of the alignment axes 
shifts. Figure 5-4 shows such a boundary where the axes exhibit a 
rotation of nearly 20“ about the vertical after correction for the 
aspect of the micrograph. 

Samples of high purity A1 were obtained and sputtered. When com- 
pared with results obtained using A1 of lesser (but uncertain) purity, 
no significant differences between the sputtered surfaces could be 
observed. Both types of Al, in the absence of deliberate seeding. 





developed undulating surfaces without a cone structure; however, small 
(several hundred Angstroms) bits of material can be seen widely scat- 
tered on the surface. These tend to be rounded and not cone-shaped. 
Possibly aluminum oxide is forming on the surface during sputtering or 
perhaps some other impurities are being deposited onto the surface. 

No cones were fomed when A1 was seeded with Ta at 475®C. 

Graphite 

Sputter time had the effect of increasing whisker density and size 
with increasing sputter time. C samples were sputtered for 5, 10 and 
20 minute intervals and showed small whisker growth, no more than 0.5 y 
in height. The etched surfaces, other than the whiskers, appeared 
slightly smoother than the non-sputtered portions of the C sample. 

Other sputter times of 1, 2, 4, and 8 hours showed a whisker length in 
approximate proportion to the sputter time. An increase in the diameter 
of individual whiskers with time was also observed along with a reduc- 
tion in the areal density of whiskers. Total sputtering time appears 
to be the most significant factor in whisker growth. 

Temperature was also a factor in whisker growth. Samples were 
sputtered for 1 and 4 hours at holder temperatures of approximately 
15°C, 300°C and 500®C. All experiments showed a decrease in whisker 
growth with higher temperatures. 

Two sets of samples were sputtered for 4 hours, one set with an 
ion beam neutralized using an immersed thermionic neutralizer and the 
other without neutralization. Various shield materials were used. All 
samples sputtered using a thermionically neutralized beam exhibited a 
greater density of whiskers that were also somewhat longer on the 


average than those on corresponding samples sputtered with a non- 
neutralized ion beam. 

yarlous shield geometries were compared: vertical edge, beveled 

edge up, beveled edge down, and circular cross section. The whiskers 
produced using a vertical shield edge or a shield with a beveled edge 
down tended to be small but uniform over the entire sample. However, 
near a shield of circular cross section or with the beveled edge up, 
the whiskers tended to cluster near the shield with a marked reduction 
in size and density further from the shield. With no shield the 
whiskers were uniformly distributed over the sample but were not 
generally as long or as dense as those closest to a beveled shield. 

Two different types of commercially available graphite were also 
compared; Union Carbide grade ATJ and POCO Graphite. No differences 
in whisker growth were observed between the two types of graphite. 

Qualitative x-ray fluorescence (an SEM attachment) and x-ray 
diffraction experiments were carried out to determine the composition 
of the C whiskers. A lattice spacing characteristic of the whiskers 
was 3.39 A. The absence of other lines from the powder pattern was a 
further indication of graphite. X-ray fluorescence had indicated Fe 
in the whisker samples. Carbon is not detectable using that apparatus 
and similar levels of Fe signal were also found using samples that 
were believed to have very little Fe. The whiskers are thus felt to 
be mostly graphite. 

Optical absorption measurements were taken of C whisker samples 
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sputtered at an etch dose of 8 mA-hrs/cm . At normal incidence and 
at a light source tilt of 64® from normal about 98.5% of the light 
from a wavelength of 3500 X to about 17000 A was absorbed. The absorp 
tlon decreased to about 97% at 25000 1 . 


Silicon 


After numerous attempts to texture Si with a variety of seed 
materials it was found that temperatures above 500®C were usually 
required. Polished, single-crystal Si wafers were sputtered at 500®C, 
550*C and 600*C for 60 minutes. Half of the wafers were cut parallel 
to the (111) plane and half were cut parallel to the (100) plane. 

Both Au and Fe were used as seed materials. 

No texture was observed on the samples sputtered using Au as the 
seed material. The following textures were thus all obtained with Fe 
seed. At 500®C the (100) Si exhibited texture only in the region up to 
about 0.5 mm from the seed source. However, at 550®C, the (100) Si 
showed dense, unifom coning up to several mm from the seed source. 

The (111) Si developed cones at 550®C and again at 600®C. The cones 
were smaller and more densely packed on the (111) Si than they were 
in the (100) Si. Beyond about 2 mm from the shield edge, the Si (100) 
600®C sample exhibited two distinct topographies, an extremely dense 
texture and a well developed cone texture (Fig. 5-5). This is the best 
texturing of Si on this scale that has been observed. The cone spacing 
averages about 1.1 p. 




VI. CONCLUDING REMARKS 

The texturing model based on surface diffusion and clustering of 
impurities has shown excellent agreement with experimental results both 
qualitatively and quantitatively. The existence of a minimum substrate 
temperature for texturing, the variation of cone distribution with 
substrate temperature and texturing with a high sputter yield seed on 
a lower sputter yield substrate were all predicted by this model. The 
model clearly has utility in facilitating the production of textured 
surfaces for a variety of applications such as low reflectance solar 
cells, heat collectors and radiators, and medical implants. The mini- 
mum critical temperature for cone formation is important in applications 
where a smooth sputtered surface is the desired end result. Because 
textured substrates are continuously cleaned by ion bombardment, the 
surface diffusion activation energies obtained should be those corres- 
ponding to a very clean surface. 

The mathematical model for sputter seeding provides a means of 
predicting and controlling the seeding rate, an important parameter in 
texturing processes. Both seeding rate and uniformity are calculable 
for a variety of experimental sputtering configurations. In addition, 
knowledge of seeding density variation across a substrate provides a 

sensitive test of the concept of a critical temperature for seed 
clustering. 

The rare gas interaction cross sections obtained will allow a 
more thorough analysis and understanding of a variety of energetic 
processes Involved in ion beam and plasma devices . Gross sections in 
this energy range (leV to iKeV) have not previously been available. 


These cross sections for elastic collision processes along with the 
previously calculated charge exchange cross sections allow designers 
of ion beam deposition and machining experiments to more readily arrive 

at appropriate configurations and pressure environments for beam pro- 
pagation. 

Experimental work on the operation of a broad beam ion source 
with CF^, which readily breaks down into reactive species, allows 
limits to be placed on the operating lifetimes of the various source 
components in a very harsh plasma environment. There are significant 
problems with thermionic electron emitters as well as with significant 
material removed from such sensitive areas as the accelerator system. 
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